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Abstract. We set up the Dirac equation in a Friedmann-Robertson- Walker geom- 
etry and separate the spatial and time variables. In the case of a closed universe, the 
spatial dependence is solved explicitly, giving rise to a discrete set of solutions. We 
compute the probability integral and analyze a space-time normalization integral. 
This analysis allows us to introduce the fermionic projector in a closed Friedmann- 
Robertson- Walker geometry and to specify its global normalization as well as its 
local form. 
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1. Introduction 

General relativity and quantum theory describe physics on different scales. Whereas 
the large-scale structure of our universe is governed by Einstein's general theory of rel- 
ativity, on the small scale gravitational forces are often negligible, and instead physical 
effects are described by quantum theory. Due to the different scales, in most situations 
it is unnecessary to combine the two theories (and in fact, a convincing "unification" of 
general relativity and quantum theory is still unknown). Nevertheless, there are issues 
in quantum theory where the large scale structure of space-time does become relevant, 
so that general relativity must be taken into account. Such problems are usually re- 
ferred to as infrared problems. To give a simple example, for the normalization of a 
wave function one must integrate the probability density over all of space, and thus 
the analysis of the normalization integral depends sensitively on the global structure 
of space. Moreover, in relativistic quantum theory, four-dimensional integrals of the 
wave functions over both space and time appear. For the evaluation of such space-time 



First author supported in part by the Deutsche Forschungsgemeinschaft. 

1 



2 



F. FINSTER AND M. REINTJES 



integrals, one must take into account the global geometry of our universe. In this pa- 
per, we shall analyze the interplay of the global geometry and integrals over quantum 
mechanical wave functions. 

Clearly, considering this problem in full generality goes beyond the scope of a single 
paper, and thus we must restrict attention to a more specific setting. First of all, we 
only consider particles of spin i, although the methods and results could easily be 
extended to scalar fields or to particles of higher spin. Furthermore, we focus on the 
physically relevant situation of the Friedmann-Robertson-Walker (FRW) geometry, 
a spatially homogeneous and isotropic space-time being the standard model for the 
present universe. In the open and flat cases the spatial volume is infinite, so that the 
normalization integrals of the wave functions will in general diverge. This is similar to 
the situation in Minkowski space, where for the proper normalization one must first 
confine the particle to finite 3-volume, for example a finite box, and can then take 
the infinite volume limit (see for example §2.6]). This construction may depend on 
the form of the finite box, whereas the global geometry of space-time does not seem 
to be relevant. For these reasons, we shall not enter such constructions here. The 
closed case is more interesting, because in this case the spatial volume and thus also 
the normalization integrals are finite. In other words, the global geometry gives rise 
to an infrared regularization. To our knowledge, this effect and its consequences have 
not yet been studied. The present paper aims at providing this analysis in reasonable 
generality and sufficient detail. 

We now introduce our problem more specifically and outline our results. In so-called 
conformal coordinates, the line element of a spatially homogeneous and isotropic space- 
time becomes 

ds 2 = S(t) 2 (dr 2 - d X 2 - f(x) 2 (d$ 2 + sin 2 dtp 2 )) . (1.1) 

Here r is a time coordinate, tp E [0, 2ir) and i? S (0, n) are angular coordinates, and x 
is a radial coordinate. In the three cases of a closed, open and flat universe, the 
function / and the range of x are given respectively by 

closed universe: f(x) = sm (x) > X £ (0,tt) 
open universe: f(x) = s i nn (x) , X > 
flat universe: f(x)=X, X > • 

The function S, the so-called scale function, is determined from the Einstein equations 
and depends on the type of matter under consideration. In the case of a perfect 
fluid, the function S must be a solution of the Friedmann equation, and the resulting 
geometry is referred to as the FRW geometry. In this paper, we do not need to specify 
the matter, and thus the scale function can be an arbitrary positive function. 
The Dirac operator on a Lorentzian manifold (M, g) reads 

V = iG j Dj , (1.2) 

where the Dirac matrices G J are related to the metric by the anti-commutation rela- 
tions 

{G j ,G k } = G j G k + G k G j = 2g jk l C 4 , (1.3) 
and Dj is the spin connection. The Dirac equation takes the form 



(V - m) * = , 



(1.4) 
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where ^ is a four-component spinor, and m is the rest mass. The spinors at every space- 
time point are endowed with an inner product of signature (2, 2) (where \& = 
is the usual adjoint spinor). Integrating this inner product over space-time, we obtain 
the bilinear form 



M 



(1.5) 



where dfiM = \f\g\d A x is the volume element on M (here g = det g^ is the determinant 
of the metric). Furthermore, choosing a space-like hypersurface TC, the probability 
integral is given by 



n 



(1.6) 



where v is the future-directed normal on Tt, and d^n is the volume element on TC. 

The main purpose of this paper is to analyze both normalization integrals (|1.5|) 
and (jl.6p for solutions of the Dirac equation in the closed universe, and to study the 
infinite volume limit S — > 00. More precisely, after a brief introduction to the Dirac 
equation in the FRW geometry and its separation (Section [2]), in Section [3] we first 
compute the spatial normalization integral (jl.6p . Then, using a WKB-type approxi- 
mation in the time dependence, we can also make sense of the space-time integral (j 1 . 5 1) 
(Proposition I3.2[) . The validity of the WKB approximation is discussed at the end of 
Section [3l These results allow us to introduce the so-called fermionic projector as a 
well-defined projection operator on the generalized negative-energy solutions of the 
Dirac equation (see Section [H (|4.2p ). In Section [5] we consider the local form of the 
fermionic projector in the infinite- volume limit and recover the Fourier integral describ- 
ing the Dirac sea in the Minkowski vacuum (see Theorem 15. ip . Finally, in Appendix El 
the spectrum and the eigenfunctions of the Dirac operator on S 3 are computed in 
detail. 



2. The Dirac Equation in the FRW Geometry and its Separation 

In this section we derive the Dirac equation in the FRW geometry and separate vari- 
ables, leaving us with an ODE describing the time dependence. The Dirac equation 
in curved space-time was first formulated by Schrodinger |15j . The more systematic 
study goes back to Brill and Wheeler (3j. Today, there are different approaches, most 
notably the formalisms using null frames [12] or pseudo-orthonormal frames pj. The 
different formalisms yield different formulas for the Dirac operator, which can be re- 
lated to each other by suitable local transformations of the spinors (see [14] for a 
discussion of this point). Here we use the approach in [5j, which is most convenient for 
concrete calculations. For an alternative derivation of the Dirac equation in the FRW 
geometry we refer to [16] . In order to satisfy the anti-commutation relations (jl.3p . we 
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choose the Dirac matrices as 
1 



QX 

G v 

where 7 , . 
sentation, 



S(r] 
1 



7 







5(r) 



(cos •& 7 3 + sin $ cos tp 7 1 + sin 1? sin 99 7 2 ) 



1 



S(r)f( X ) 
1 



sin $ 7 3 + cos •& cos tp 7 1 + cos # sin </? 7 2 ) 



sin 99 7 1 + cos (p 7 2 ) 



(2.1) 



S(r)f(x) sintf 

, 7 3 are the usual Dirac matrices of Minkowski space in the Dirac repre- 



7 





-1 



7 



and a a , a = 1,2, 3, are the Pauli matrices, 



a 



1 

1 



a' 



-i 

1 







a 



a 



The spin connection Dj can be written as 



Dj = 8a — iE-j 



(2.2) 



(2.3) 



Sj — Wj " J J 

where Ej, the so-called spin coefficients, are given in terms of the Dirac matrices and 
their first derivatives. Substituting f|2.3|) into (jl.2|) . one sees that in order to obtain the 
Dirac operator, it suffices to compute the combination G^Ej. In the FRW geometry, 
this combination simplifies considerably: 

Lemma 2.1. For the choice of the Dirac matrices 112.1]) , the spin coefficients satisfy 
the relation 

&E 3 = % -V 3 G\ (2.4) 
where V denotes the Levi- Civita- connection. 
Proof. According to [5], the spin coefficients can be written as 



E; 



p{d jP ) 



4 Tr(G m V 3 G n )G m G n + 1 Ti{ P G 3 V m G m ) p 
lb o 



(2.5) 



where 



P ■-- 



i 

4! 



\g\ tijki G l G J G k G l , 



and eijki is the totally antisymmetric symbol with £ T x$<p = 1- 

Let us show that the first summand in (|2.5I) vanishes. Since the FRW metric is 
diagonal, different Dirac matrices anti-commute, i.e. G l G 3 = —G 3 G l if i 7^ j. Using 
furthermore the anti-symmetry of the e-symbol, we conclude that 

p = i x /\^\G T G x G®G v . (2.6) 

Substituting the formula \/\g\ = S' 4 (r)/ 2 (x) sin^? and using the definitions (12. lh . a 
straightforward calculation shows that p = Z7°7 1 7 2 7 3 . Hence p is a constant, and so 
its partial derivative in the first term in (12. 5ft vanishes. 

Next we show that the last summand in (|2.5p vanishes. Since the G- 7 are linear 
combinations of the , i.e. G n = a n j r y J with real coefficients a™, we find 

V n G n = V„(a n i ) 7 j , (2.7) 
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and thus 

TripGjVmG™) = aft m (a m k ) Tr(p 7a k ) ■ 

Due to the anti-commutation relations, the last trace vanishes for any choice of Dirac 
matrices 7* and 7*, proving that the last summand in (|2.5p indeed vanishes. 
Substituting the remaining second summand in (|2.5p into (|2.4p . we obtain 



16 



Tr(G m VjG n ) G j G m G n , 



(2- 



where we raised and lowered the indices with the metric. In the case where the 
indices m, j, and n are all different, the expression is totally anti-symmetric in these 
indices, because the last three Dirac matrices in (|2.8p anti-commute. Consequently, in 
this case the symmetries of the Christoffel symbols allow us to replace the covariant 
derivative by a partial derivative, and an explicit computation of the partial derivatives 
shows that the trace vanishes. Hence it remains to consider the case that at least two 
of the indices m, j, and n coincide. Ricci's lemma and the Leibniz rule yield that 

= 4V jg mn = VjTr(G m G n ) = Tr(V jG m G n ) + Ti(G m VjG n ) , 



and thus we may anti-symmetrize the expression (|2.8p in the indices j and m. This 
gives rise to the simplification 



GPEi 



TT{G n V m G m )G n . 



According to (|2.7p . the expression V m G m can be written as a linear combination of 
the Dirac matrices 7 J , or alternatively of the matrices G 3 . Applying the relation 



Tr(G n G j )G n = 4g nj G n = 4G j 



gives the result. 



□ 



A major advantage of the relation (|2.4p is that it involves a divergence, which can 
be computed using the Koszul formula 

1 



V m G" 



\9\ 



g\ G r ' 



making it unnecessary to compute the spin connection coefficients or even the Christof- 
fel symbols. A short calculation combining Lemma 12 . 1 1 with equations (|2,3p and (jl.2p 
yields for the Dirac equation in the FRW geometry 



iG T d T + 



35 
2S 



+ iG x [d x + 



f'-l 
f 



+ iG»d# + iG^d^ 



m 



* = 0, 



(2.9) 



where the dot and prime denote the partial derivatives with respect to r and Xi re_ 
spectively. Multiplying by S(t), this equation can be written more conveniently as 



S(r) m + 



V n 
-Vn 



* = 



where the purely spatial operator V7-1 is given by 



d Y + 



f'-i 



(2.10) 



(2.11) 
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and the matrices a a , a G {x,$,*p}, are linear combinations of the Pauli matrices, 
* x •= cos i9 <t 3 + sin $ cos ip a 1 + sin $ sin ip a 2 







1 



(7 



fix) 



a 



<fi .- 



1 



sin ■& a 3 + cos ■& cos <p a 1 + cos # sin <r 2 



sin (p a 1 + cos 93 <r 2 ) . 



(2.12) 



/(x) sin?? 

In the case of a closed, open and flat universe, the operator D-^ is the Dirac operator 
on the sphere S 3 , on an hyperboloid, and on flat M 3 , respectively. From now on, we 
shall restrict attention to the closed case, which has the advantage that the spectrum 
of T>-ft is discrete. As is worked out in detail in Appendix [Al the operator can be 
identified with the intrinsic Dirac operator on S 3 . It is an essentially self-adjoint elliptic 
operator on the Hilbert space L 2 (S 3 ) 2 with domain of definition C°°(S 3 ) 2 (see |llj). 
It has the purely discrete spectrum [11] 



and the dimension of the corresponding eigenspaces is 



± 3 - ± 5 - ± 7 - 

2' 2' 2' 



(2.13) 



dimker(P w - A) = A 2 



1 

4 



In Appendix|A]an orthonormal eigenvector basis is given explicitly in terms of spherical 
harmonics and Jacobi polynomials (see Theorem IA.5|) . It is denoted by (Vi-t)j where 



n G Nq, j € Nq + \ and k G {—j, + The eigenvalues are given by 



with 



A = ±(n + j + l) 



(2.14) 



njk ^njk 

More generally, we denote a normalized eigenfunction of T>~n corresponding to the 
eigenvalue A by tp\ G L 2 (5 3 ) 2 . 

Employing for ^ the separation ansatz 



S(t) 



(2.15) 



we obtain a coupled system of ODEs for the complex- valued functions hi and hi, 



id T — Sm 







+ A 



1 

1 



. 



(2.16) 



3. Normalization Integrals in the Closed FRW Geometry 



In this section we shall analyze the integrals (|1.5j) and (|1.6p for our separated wave 
function (j2. 15f) . To compute the probability integral (|1.6I) . we choose TL to be a slice 
of constant conformal time r. Then the future-directed normal v has the components 
v T = S(t) and v a = for all a G {x,$,p}- The volume element on TL is dfin = 
\/\gn\ d>x d"d dip = S 3 (r)dn S 3, where is the volume element on the unit sphere S 3 
(thus g-j-i is the determinant of the induced metric) . Hence 



\E ,t \I> dji 



n 



H 



I 2 + N 2 ) IV'aI 2 ^ 



(3.1) 



From the ODE (12. 16ft one sees (using that the matrices in (12.16P are all Hermitian) that 
the function |/ii| 2 + |/i2| 2 is constant in time. This corresponds to the more general fact 
that, as a consequence of current conservation, the probability integral is independent 
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of the choice of TC. Adopting the same convention as in [61 §2.6], we normalize the 
wave functions such that 



(3.2) 



The volume element in the space-time integral (|1 .5|) is given by 
dfJ-M = \f\g\ dr dx d'd dtp = Sdr d/j,-n . 
Thus for two wave functions ^ and fy, 

<^|\J/>= / dfi M = / Sdr / dn n ~$ii! 
Jm J Jn 



Sdr J dfi H (hihi - h 2 h^ (ip x ,^ x )c 2 > 
where the bar denotes complex conjugation. Hence 

<*|*> = W>a,-0a>l 2 (53)2 / (Mi-^Aj) Sdr . (3.3) 



We first note that if the universe has a finite life time, then this normalization 
integral is necessarily finite. 

Proposition 3.1. Consider a spatially homogeneous and isotropic geometry (jl.ip for 
the conformal time in the range T m \ n < r < r max and singularities at r m \ n (the "big 
bang") and at r max (the "big crunch"). Normalizing the wave functions according 
to (|3,2p . the space-time inner product (jl.5p is finite and bounded by 

|<tf I tf>| < — / Sdr. 

° 'mm 

Proof. Working in (|3.3[) with normalized spatial eigenfunctions, we obtain 
<tf|tf></ \\hi\ 2 - \h 2 \ 2 \ Sdr < / + \h 2 \ 2 ) SdT . 

^ Tmin ** ^"min 



The result follows immediately from (|3.ip and the normalization convention (|3.2p . □ 



Clearly, the time integral in (|3.3p is very large and will in general diverge in the 
limit when the life time of the universe tends to infinity. In [6, §2.6] this problem is 
bypassed by working with a variable mass parameter. In order to get the connection 
to this normalization method, for a fixed spatial eigenvector ip x we now consider so- 
lutions ty m and \E' m of the form (|2.15p for two variable mass parameters m, m! > 0. 
Then the corresponding time-dependent functions (/i™,/i™) and (h™' , h™' ) are solu- 
tions of (|2,10p . We want to compute the inner product <*$> m \ty m ' > , which according 
to (13.31) becomes 



<i$ ,m^m' > = {Whf -K?hf'\SdT . (3.4) 

" Tmin 

Since the Compton wave length is much smaller than the life time of our universe, the 
integrand will typically be highly oscillatory. Therefore, it seems appropriate to use 
a WKB-type approximation (this approximation and its limitations will be discussed 
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in detail at the end of this section). To this end, we diagonalize the matrix potential 
in (|2.16p with a unitary matrix U(t), 



U 



Sm 
-A 



-A 
-Sm 



) U- 1 = Vm 2 S 2 + A 2 (I 



Then (|2.16p can be written as 



id t 



U 



Vm 2 S 2 + A 2 



1 
-1 



U 



(3.5) 



(3.6) 



Under realistic conditions, the term U is very small compared to mS. Leaving out this 
term, we can solve the ODE explicitly by 



(&)m-"<t)-' 



/ 



ci exp J \J m 2 S 2 + A 2 dr^j ^ 
\J m 2 S 2 + A 2 dr 



c 2 exp 



(3.7) 



This approximation has the nice property that it respects current conservation \hi\ 2 + 
I h-2 1 2 = const, and thus the normalization f|3.2|) is implemented simply by the condition 



I 1 2 I |2 
Fil +l c 2| 



1 

2t7 



(3.8) 



Indeed, this property is why we prefer (13. 7j) over the alternative method of rewriting the 
Dirac equation in terms of scalar second order equations and employing the standard 
WKB ansatz. 

Substituting the WKB ansatz (|3.7p into (|3.4h and multiplying out, we obtain inte- 
grals of the form 



/ 77(7-) exp ji/ (=F \Jm 2 S 2 + A 2 ± V m' 2 S 2 + X 2 ) df 1 dr , 

Jr min v J To > 



(3.9) 



where tq determines the phase of the exponential, and 77 stands for the integration 
density S times a combination of the constants c^L , and certain matrix elements of 

the unitary transformations U m and U m . Qualitatively speaking, the functions 77(7") 
and S(t) vary only on the cosmological scale, whereas the exponential oscillates on the 
microscopic scale, unless the curly brackets in (|3.9|) vanish. As a consequence of these 
rapid oscillations, the value of the integral will be small, except if the square roots 
terms in (|3.9|) have opposite signs and m ~ m! . In order to quantify the contribution 
in this limit, in (|3.9f) we only take the linear term in 8m := m — m! to obtain 

I := I ^ t](t) exp U I U/m 2 S 2 + \ 2 - Vm' 2 S 2 + A 2 ) df\ dr 

J T min I J TQ ) 



77 (r) exp < i 5m 



mS 2 



TO 



Vm 2 S 2 + A 2 



df + 0((<5m) 2 ) \dT 



(and similarly for the other combinations of signs). After multiplying by a convergence 

for any fixed r) over 5m using the forn 

where Q = [ , df . 



generating factor e £ ( 5m ) 2 ; we can integrate (for any fixed r) over 5m using the formula 

mS 2 



JQSm-e(6rn)i 



7T _Ql 

e 4 e 



Vm 2 S 2 + A 2 
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We thus obtain to leading order in 8m, 

f I e~ £ ^ 2 d(6m) = / TmaX n(r) exp dr . (3.10) 

In the limit e \ 0, the exponential tends to zero unless Q vanishes, and thus we may 
evaluate the integral in the saddle point approximation. Since the only zero of Q is 
at r = To, we may replace Q in (|3.10p by 

Q « Q'(ro) (r - T ) , Q'(to) = . = 2 

y / m 2 S(T ) 2 + A 2 

to obtain 

I e-<^ 2 d(8m) = P" V(r ) ^ dr ^ 2vr . 

V e A min <2'Oo) 

We conclude that in the above saddle point approximation, the oscillatory integral (|3.9p 
can be computed by 

/ ??(T)exp<i / (sy/ m 2 S 2 + A 2 — s'y m' 2 S 2 + A 2 J df > dr 
Jr min I J to ) 



2vr 6 StS , 5(m - m') r](r ) V , (3.11) 



V m 2 5(T ) 2 + A 2 
m5(r ) 

where s,s' £ {1, —1}. We remark that this result could be obtained more directly if 
one replaced the functions r\ and S in (|3.9p by the constants rj{jo) and S'(to), replaced 
the integral from T m ; n to T max by an integral over the whole real line and applied the 
well-known distributional formula 

Au)t 



e luJt = 2tt 8(u>) . 

However, our method clarifies that the phase of the exponential in (13, 9p as described 
by the parameter to determines at what time the functions r\ and S in (|3.1ip are to 
be evaluated. 

Evaluating the integral (|3.4p with the help of (|3.1ip gives the following result. 

Proposition 3.2. For a family of wave functions (^ m ) m >o °f the form (|2.15p with a 
WKB ansatz for the time dependence (|3.7p . the inner product (jl.5p can be computed 
in the saddle point approximation to be 

<q> m | = 27r £( m _ m /j (| Ci |2 _ | C2 |2^ _ (3 12 ) 

Proof Substituting (|3.7j) into (|3.4p . according to (|3.1ip the mixed terms ~ c™c™' 
and ~ d^c™' do not contribute. Thus 

r T max f j \ 

<q ,m \^m >= / | c m c m ^m^r/™ )\ - (U m ) 2 (U m )?J e"** 

■* ^~min 

+ W (W™)l(U m ')l - W%{U m ')t) e**} S'dT, 
where we used the abbreviation 

$= / (Vm 2 S 2 + A 2 - V™' 2 ^ 2 + A 2 ) df . 
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Again using (|3.1ip . we only get a contribution at m = m', 

2vr5(m - to') + A ' 

to6 

ci| 2 W\ 2 - |tf 2 | 2 ) - |c 2 | 2 (|C/f | 2 - |C#| 2 ) ), (3.13) 

where we omitted the arguments To and the indices m. The unitary transformation U 
in (13. 5h can be given explictly by 

C/ = cosg)l + sing)(j - 1 ) with r = arctan(A). 



Thus 



|[/1|2 _ |[/2|2 = |^2|2 _ |^1|2 = cos 2 ' \ _ sin 2 , _ ) = cog T 



mS 



■2J \2J Vm 2 S 2 + A 2 ' 

giving the result. □ 

We conclude this section by explaining the result of Proposition 13.21 and discussing 
its physical significance. Considering a variable mass parameter together with the 
5{m — m')-normalization (I3,12|) is very useful for the construction of wave functions 
for which the space-time integral (|1.5p is finite. More precisely, choosing ij G C^°((m — 
e, m + e)), we can "smear out" the mass parameter on the scale e by setting 



^ = j n(m) ^ m dm (3.14) 

(more generally, one can apply this method to composite expressions, see for exam- 
ple flO}). Then 



<^r|^r>= dmr](m) I dm' r](m') <^ m > = \r](m)\ 2 dm . 



The surprising point about Proposition 13.21 is that the scale function S(t) drops 
out when computing (|3.13p . making (|3.12p independent of the space-time geometry. 
In particular, the parameter tq describing the mass dependence of the phase in (|3.9p 
does not enter the 5{m — m')-normalization. In this way, the result of Proposition 13.2 
is robust to the geometry of our homogeneous isotropic space-time and to the phases 
of the function r\ (to) in the mass smearing (|3,14p . 

Physically, the mass smearing (|3.14|) can be understood as a technical device for 
introducing a finite length scale ~ e _1 for the correlation of \& with itself. Such a finite 
correlation length could also arise as a consequence of an interaction or self-interaction 
of a fully quantized system. In view of this more general picture, the conclusion of the 
above construction is that if the correlation length is made finite in a Lorentz invariant 
way, then the space-time integral (jl.5p should be controlled explicitly in terms of the 
probability integral (11. 6p . 

Finally, the WKB approximation (13. 7p requires a detailed explanation. To avoid 
misunderstandings, we first point out that we used the WKB ansatz only in the time 
dependence, whereas the spatial dependence was treated without any approximations. 
Consequently, the WKB approximation only enters the analysis of the time integral in 
the space-time inner product (|1.5j) . giving the result of Proposition l331 The probability 
integral (|1.6p . however, was computed exactly (see (13.21) ). The time independence 
of the probability integral is a general consequence of current conservation, without 
referring to any approximations. 
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Qualitatively speaking, the WKB wave function is a good approximation if the 
potential is almost constant and non-zero (see for example [9]). More precisely, com- 
puting the error term in (|3.6p . one finds that the ansatz (|3.7p is a good approximation 
provided that the following condition holds: 

|A|m 3 S«1. (3.15) 
(m 2 5 2 + A 2 )5 

For our present universe, the expression on the left is extremely small. Namely, we 
can estimate it by 

|A|m 5<4,-2-10- 38 



(m 2 S 2 + X 2 )l mS2 

where we used that the quantity S/S 2 is the Hubble constant H « 2, 3 ■ 10 s , 
whereas for m we set one over the Compton wave length of the electron, m ~ 1.24 • 
10 20 s -1 . We conclude that at present, (|3.7h is an extremely good approximation 
describing fermions which move freely in the FRW geometry. Clearly, the WKB ansatz 
does not take into account particular quantum effects like pair creation (because the 
components of positive and negative frequency are not "mixed" in the WKB ansatz; 
for a more detailed discussion of this point we refer to [TO]). But since we here focus 
on the gravitational interaction, such quantum effects will be very small and should 
not be of relevance for all normalization issues. Furthermore, the WKB approximation 
does not allow for the description of particular quantum effects in the gravitational 
field like the Unruh or Hawking radiation (see for example [17j). However, these effects 
do not appear in the FRW geometry and can thus be disregarded. 

Another shortcoming of the WKB ansatz (|3.7p is that it does not apply near the 
big bang or big crunch singularities. Indeed, as shown in [8], near these space-time 
singularities the Dirac spinors show quantum oscillations, which cannot be described 
semi-classically and lead to unexpected quantum effects which can even prevent the 
formation of the singularity. Even in such rather wild quantum scenarios, the WKB 
approximation is of use and allows us to compute the space-time integral (jl.5p . as we 
now explain. Proposition 13.11 (which does not use any approximations) tells us that 
the space-time integral is a-priori bounded by J Sdr. Thus we know that the time 
intervals shortly after the big crunch and shortly before the big bang yield only very 
small contributions to the inner product (]1.5p (no matter how wild the behavior of the 
Dirac spinors and the metric is). Away from the singularities, we can again use the 
WKB approximation to obtain the result of Proposition 13.21 

We conclude that, although the WKB approximation does certainly not account 
for all physical effects of Dirac spinors in curved space-time, for the purpose of ana- 
lyzing the normalization integrals on the large scale, the WKB ansatz (|3.7p seems an 
admissible and indeed very good approximation. 

4. The Global Normalization of the Fermionic Projector 

We now turn attention to the normalization of the Dirac sea in the FRW geometry. 
In the vacuum, the Dirac sea is composed of all negative-energy solutions of the Dirac 
equation. In our time-dependent setting, the energy of the Dirac particles is not 
preserved in time, and thus it is not obvious of which states the Dirac sea should 
be built up. In [H [7] this problem is resolved for a general time-dependent external 
field by a global construction in space-time, giving a canonical splitting of the solution 
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space into two subspaces, one of which generalizes the negative-energy solutions of 
the vacuum and thus extends the concept of the Dirac sea to the interacting situation. 
However, because of technical assumptions on the decay of the external field at infinity, 
these results do not immediately apply in the FRW metric, and the constructions 
of [H [7] have not yet been extended to the FRW geometry. Fortunately, the problem 
of distinguishing the "generalized negative energy solutions" disappears in the WKB 
approximation (|3.7p . because then the negative-energy part is obtained simply by 
setting ci = 0. Working again with the approximation (|3.7|) . we can describe the 
Dirac sea explicitly, as we now explain. 

Following [U [7J , the fermionic projector P is an operator on the Dirac wave functions 
whose image is spanned by the states of the Dirac sea. Furthermore, the fermionic 
projector should be symmetric and idempotent with respect to the inner product (jl.5p . 
Thus, using a bra/ket-notation, 

P = -El^><*A, ff |, (4.1) 

A.<7 

where A refers to the eigenvalues of the spatial Dirac operator P53, and a labels a basis 
of the corresponding eigenspaces of T> S 3, and the Dirac wave functions must satisfy 
the normalization condition 

It is convenient to write (|4.1j) with an integral kernel, 



(P*)(x) = / P(x, y) V(y) duM with P(x, y) = - V %xA x ) *A,«r(l/) • 

Employing the ansatz (|2.15p and working for the time dependence with the WKB 
approximation (|3.7p . we select the generalized negative-energy solutions by choos- 
ing ci = 0. In order to satisfy (|3.8p . we choose C2 = (2ir)~2 . Clarifying the dependence 
on the mass parameter by an index m, we thus obtain for the kernel of the fermionic 
projector 

A,cr 



?) U (y°) (J ex P (*£ Vm^S^ + X^dT^j , (4.2) 



where the ip\ j(7 form an orthonormal basis of the A-eigenspace of P53 and x°, y° denote 
the time coordinates. 

Clearly, the above operator P m is symmetric with respect to the inner product (|1.5p . 
Furthermore, using the normalization results (|3.8p and (|3.12p . it is idempotent if we 
work with a ^-normalization, 

P m P m > = $( m ~ m ') p m ■ 

This ^-normalization can be avoided similar to (13.140 by "smearing out" the mass 
parameter, for example by setting for any 5 > 

cm+S 



pm+o 
P= P^djl 
J m—S 
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Then this P really is idempotent, 

rm+5 rm+8 rm+5 

P 2 = dn\ dfx' P m P m , = P m = P. (4.3) 

J m— S Jm—5 Jm—5 

5. The Local Form of the Fermionic Projector 



We shall now explore how the fermionic projector (|4.2p looks like locally, in the 
physically realistic case when the size of the universe is much larger than the length 
scale for local observations. Our analysis involves the two approximations that the 
scale function S is very large and that it is almost constant locally. 

Theorem 5.1. We consider the fermionic projector (|4.ip with the individual states 
normalized according to (|3.2p . Then in the limit S — > oo of a large universe and 
neglecting local curvature effects, the fermionic projector (|4.ip in a local reference 
frame takes the standard form of Minkowski space, 

Pm(x, y)= I t~~tt {H + m) 5(k 2 - m 2 ) @{-k°) e'^'^ . (5.1) 



(2ny 

Proof. For convenience, we choose the local reference frame such that y = 0. In order 
to bring the Fourier integral in (|5,ip into a more convenient form, we set to = k°, 
t = x , r = |x| and choose polar coordinates (p = |A;|,i?, <p) in momentum space with 
the north pole pointing in the direction of x (thus i? denotes the angle between k 
and x). We then obtain 



/ 



dk (# + m) 5(k 2 - m 2 ) e{-k°) e~ lkx 



(2vr) 4 



1 f'OG poo rl pZ7T 

— ^i(i@ x + m) / duj I p 2 dp dcosti / dip 5(u; 2 - p 2 - m 2 ) e~ iujt+iprcos{> 
( 27r ) J-oo Jo J-i Jo 



1 , f°° P 2 dp rl 



(2tt)3 ^ ;j o 2 \ u(p)] 



(i@ x + m) / rfy^e-^ (p)t / e iprcos ^ dcostf 



i 



1 HO,. m) [ X JUllL^l, i 



(2tt)3 ^* >J |u,(p)| 

1 f°° p ( , , . o . sin(pr) . / sin(pr) 



(2tt)3 J " | W (p)| 



* HPh° + m) + vfDr e-*«* , (5.2) 



where we set oj(p) = —\/p 2 + m 2 and 



7 r = cos i? 7 3 + sin i? cos ip 7 1 + sin t? sin 99 7 2 . (5-3) 

Our goal is reproduce (|5.2p from ()4.2p by suitable approximations. First, neglecting 
local curvature effects, we may assume that S is constant in our neighborhood. Fur- 
thermore, we may disregard the r-dependence of the transformation U. From (|3.5p 
one sees that 

v -x ( °) u _ 1 1 fSm " A 



V " 2 + A 2 V" A ~ Sm 
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Moreover, rewriting the spatial eigenfunctions in (|4.2p in terms of the spectral projec- 
tors E\ as computed in the appendix (see (|A.24|) and (|A.25|) ). we obtain 

P m (x,y) 



1 



1 



2ttS 3 



where 



to 







A 

+ S 





-1 



+ m 



s iwS(y°-x°) 



(5.4) 



A2 
S 2 



+ m z 



In order to further simplify this formula, we fix y and choose our coordinate system 
such that y° = and that the point y E S 3 is at the north pole. Describing the 
point x by the coordinates (r, x> <p), the spatial spectral projectors can be expressed 
by Lemma lA.61 Next, it is important to observe that the quantities r, \ an d A are 
dimensionless. In order to obtain the usual physical observables time, radius and 
momentum, we must multiply by suitable powers of S, 



t = St 



Ssinx » Sx , 



|A| 



These local observables should be well-defined in the limit S — > oo. This means 
in particular that the spatial eigenvalues |A| must tend to infinity. This allows us 
to simplify the formula of Lemma IA.6I using the large-n-asymptotics of the Jacobi 
polynomials [TJ eq. 22.15.1] and the Stirling formula [TJ eq. 6.1.37]. This gives 



E ±in+ 3 ) (x,y) = ^ 



nsin(nx) 



ncos(nx) sin(nx) 



I +C)| - 

n 



X V X X 

Substituting this formula in (15. 4p . and multiplying out, we can rewrite the tensor 
products of 2 x 2-matrices in terms of Dirac matrices (also compare (I2.12P and (|5.3p ). 



1 



C 2 





-1 



7 





-1 



7 



Moreover, inspecting the signs in the case of positive and negative A, one sees that 
both cases give the same contributions. Thus the fermionic projector becomes 



P m (x,y) 



1 



8vr 3 5 3 



En 
\r,A 



raeNo 



, xSin(nx) . r 
+ m) h ipy 

X 



cos(nx) sin(nx) 



X 



nx 



up to corrections of higher order in 1/n. Using that |A| = n + I, we can set p = n/S. 
Furthermore, considering the n-series as a Riemann sum, one sees that the sum goes 
over to an integral via 

1 \ - S^oo 



nGN 



dp ■ 



(5.5) 



and at the same time the corrections of higher order in 1/n tend to zero. We thus 
recover (15.21). □ 



We point out that the same normalization constant was obtained in [U §2.6] by 
considering the system in finite 3-volume in Minkowski space and taking the infinite- 
volume limit. We here obtain the same result using a more realistic infrared regulariza- 
tion in the form of the closed FRW geometry. Our result shows that despite the naive 
divergence of the space-time integral (|1.5p . the fermionic projector is well-defined in 
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the infinite volume limit. Our analysis also confirms that the normalization constant 
is independent of the details of the regularization procedure. 

We remark that the above methods could also be used to analyze in detail how 
the fermionic projector in the FRW geometry deviates from that of Minkowski space. 
Namely, the local curvature is captured by the Jacobi polynomials. The effect of the 
global geometry could be analyzed by considering the difference of the integral and 
the Riemann sum in (|5.5p . However, this analysis would go beyond the scope of this 
paper. 

6. Conclusion 

We saw that the Dirac equation in the closed FRW geometry can be separated 
into ODEs describing the spatial and the time dependence. The spatial dependence 
can be solved in closed form (see Appendix [A]). This makes it possible to compute 
the probability integral (jl.6p . which is time independent due to current conservation. 
For the space-time normalization integral (|1.5p the situation is more difficult, because 
the solution of the time-dependent ODE cannot be given in closed form. Thus we 
must rely on estimates (Lemma 13. If) as well as on the WKB approximation (13, 7p . 
which are justified at the end of Sectional After "smearing out" the mass parameter 
(see (|3.14p ) , we can make sense of the space-time integral (ll.5p (Proposition I3.2f) . We 
apply these results to construct the so-called fermionic projector, being a well-defined 
projection operator on the generalized negative-energy solutions of the Dirac equation 
(see (14.2p ). Locally, the fermionic projector has the same form as in Minkowski space 
(see Theorem I5.ip . In this way, a proper normalization of the fermionic states was 
obtained in a physically realistic setting. We also learn that the details of the global 
geometry do affect neither the local form of the fermionic projector nor the relation 
between the normalization integrals (II. 5p and (II. 6p . 

Appendix A. Spectrum and Eigenfunctions of the Dirac Operator on S 3 

In this appendix we derive the intrinsic Dirac operator on S 3 (with the standard 
metric), compute its eigenvalues and construct an explicit orthonormal eigenvector 
basis in terms of special functions. We also derive convenient formulas for the spectral 
projectors. Clearly, the results of this section are not new. But it seems worth to give 
a self-contained and explicit treatment, in particular because the detailed formulas 
for the spectral projectors are needed in Section [5j As in (jl.ip we choose coordi- 
nates (x, £ (O; 71 ") x (0,7r) x (0, 2-7r), such that the line element on S 3 becomes 

ds 2 = dx 2 + sin 2 x d'd 2 + sin 2 x sin 2 i? dip 2 . 

We first verify that the operator (|2.1ip really is the intrinsic Dirac operator. 

Lemma A.l. The intrinsic Dirac operator on S 3 can be written as 

V s3 = ia x (d x + COSX ~ 1 ^ + ia^dv + y<>d v , (A.l) 
where the matrices a x , o~ and are given by \2.1S\) . 

Proof. Before we can apply the methods of [5], we must get back to a Lorentzian 
manifold. To this end, we consider the manifold M = Rx S 3 with the line element ds 2 
given by 

ds 2 = dr 2 - ds 2 , (A.2) 
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where r £ M is the "time coordinate" and ds 2 again denotes the metric on 5 3 . Then 
the hypersurface r = is isometric to S 3 , and the corresponding second fundamental 
form vanishes identically. As a consequence, the Dirac operator on M can be written 
as a constant time derivative term plus a term involving the intrinsic Dirac operator 
on S 3 , 

The line element (1A.2P is obtained from that of the closed FRW metric (II. ip by set- 
ting S = 1. In this special case, the Dirac operator in (|2.9p reduces to the Dirac 
operator in (I2,10p with S = 1. Comparing with (IA.3D gives the result. □ 

Our first step in diagonalizing the operator (jA.ip is to separate the and in- 
dependence in (jA.ip . The next lemma gives a connection to the standard angular 
momentum operators 

L = — ix A V and K = aL + 1 , (A.4) 
where "A" denotes the cross product in R 3 . 

Lemma A. 2. The angular Dirac matrices in (|2.12p and the angular momentum op- 
erators (|A.4p satisfy the relations 

a°d# + a v d v = aL (A.5) 

Ka x = _ a x K _ ( A 6 ) 

Proof. Choosing in Euclidean M 3 the polar coordinates (sinx, i9, y?) with sin^ = \x\ , 
a short calculation shows that 



o-®d$ + a^du, = cxV — a x d s i nx = — — ax (aV — a x d S i ny 

sinx v 

(ax)(aV) = xV + ia{x A V) = sinx<9 s i nx — aL . 

Combining these relations gives (|A.5[) . Furthermore, using the Leibniz rule as well as 
the anti-commutation relations of the Pauli matrices, we obtain 

sinx (aLa x + a x aL^j = -i jc?(x A V), rz^cr" j = -2ix(xf\\7)-ia(xf\a) = -ia(xAa), 

where we used the Einstein summation convention for the index a. Writing the cross 
product with the totally antisymmetric Levi-civita symbol e Q ^ 7 and using the relation 

a a a 13 = ie^a 1 , 

we find 

-ia{x Aa) = -i e a ^a^x a a fi = x a e a ^e^ s a s = -2xa , 

and thus 

aL a x + a x aL = . 

smx 

Using the definition of K in (|A.4p we obtain (|A.6P . □ 

Following the procedure in quantum mechanical textbooks, we next diagonalize the 
operator K in terms of the spherical harmonics . 
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Lemma A. 3. The two- component wave functions 6 L 2 (S 2 ) 2 with j 6 No + \ 

and k £ {-j, -j + 1, j} given by 



X JH \ " 2, ■ 2 Y i+i [o) f 2j + 2 VI 



i + 1 - k fc-i /A H + 1 + k fc+i /0 



(A.7) 



/orm an orthonormal eigenvector basis of the operator K in (|A.4p . urai/i the eigenvalues 
given by 

^k^MHr (A - 8) 



Furthermore, 

"2 J^2 



o- x xLl = X-±l ■ (A.9) 



Proof. Following standard conventions, the spherical harmonics YM'd, if) G L 2 (S 2 ) 
with / 6 No and £ {-1,-1 + 1,...,!} are an eigenvector basis of the angular momen- 
tum operators L 2 and L 2 , 

l 2 ^ = /(/ + i)y/= , L z y^ fe = krf . 

As a consequence, the "ladder operators" L± := L x ± iL y satisfy the relations 

L±Y? = + 1) - k(k ± 1) Y/^ 1 , 

where for convenience we used the convention = if |/c| > /. Using these relations 
together with the explicit form of the Pauli matrices (|2.2j) , a straightforward calculation 
gives (|A.8|) (for more details see [13J). The orthonormality of the wave functions (|A.7|) 
is verified directly using that the spherical harmonics are orthonormal in L 2 (S 2 ). The 
completeness of the x k +i follows because taking suitable linear combinations of the 

^ ± 2 

wave functions in (|A.7|) we obtain the spinors 

y*7(T) and y^'f? 

which clearly form a basis of L 2 (S 2 ) 2 . 

In order to explain (jA9]), we first apply (|A.6|) to obtain 

J^y* i = -or x K X k i = T ( i + - ) ^ X X^ i • 

Hence the operator a x maps the eigenspaces of K corresponding to the eigenval- 
ues ±(j + i) into each other. Since this operator is unitary, it is clear that the 
vectors (o~ x x k i)k=—j...j form an orthonormal basis of the eigenspace corresponding 

to the eigenvalue =f(j + ^). This proves (|A.9|) up to unitary transformations of the 
eigenspaces. To verify that these unitary transformations are the identity, one needs 
to go through a straightforward computation using the detailed form of the wave func- 
tions (IA.7j) and of the matrix a x as given in (|2.12p , □ 

We are now ready to separate the Dirac equation on S 13 . In the Dirac equation 

V S 3ip = Xtp 
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we substitute (jA.ip and multiply from the left by —ia x . Using (|A.5|) and (|A.4p . the 
Dirac equation becomes 

V sinx smx / 

Next, for any j € No + \ and k G {—j, + we take the ansatz 

V(x, &,<p) = ^- (*i(x) X i " »*2(X) xj+ito ¥>)) (A.10) 

with two complex functions $1 and $2- Applying (IA.8jl and (jA.9[) . we obtain the 
system of ODEs for the two-spinor $ = ($1, $ 2 )> 



1 



8x _^±l 1 +A 1 
x smx V° _1 7 V _1 



$ = 0. (A.ll) 



We refer to this system as the radial equations. It can also be written as the eigenvalue 
problem 

7» = A* where J) 4, + £+1 (° J) . (A.12) 

Using the orthonormality of the eigenfunctions x k > 1 > we find that 

^,^) i2(5 3 )2 = / \$>\ 2 d X - (A.13) 

JO 

Hence our goal is to find all solutions to (|A.12p which are normalizable with respect 
to the standard L 2 scalar product on the interval (0, 7r). 

In the special case A = 0, the ODEs (jA.lip can easily be solved in closed form to 
obtain the two fundamental solutions 

Every non-trivial linear combination of these two function has a non-square-integrable 
singularity at x = or X = 71 j an d thus A = is not an eigenvalue. In the remaining 
case A 7^ 0, we can solve the second equation in (jA.lip for $1, 

*i = T ( $ 2 + — *2 ] • (A.14) 
A y smx J 

Differentiating this relation and substituting both 3>i and in the first equation 
in (|A.lip . we obtain the second order scalar equation 



+ ' ( C ° S X + j + \ I + A ' 1 $2 = ° ' (A,15) 



We next construct the general solution to the ODEs (|A. 15|) . We take the ansatz 
$2(x) = Kx)g(x) with 



1 — cos x 

and perform the coordinate transformation 



1 + cos x\ 2 



x^y-= — ^ — • (A.16) 
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This leads to the hypergeometric differential equation 

(y 2 - y) g"(y) + {y + (2j + 1)) g'(y) - A 2 g(y) = o , 

whose general solution can be written in terms of the hypergeometric function 2 Pi by 

g(y) = a 2 F 1 (-\,X,-j,y) + (3 y^ +1 2 F 1 (-\ + j + 1, A + j + 1, 2 + j, y) (A.17) 

with two complex parameters a and (5. We conclude that the general solution to (|A. 15|) 
is 

s i+i , 
1 + cos x \ 2 ( 1 — cos x 



•»<*>- ^-t-m- (A - 18) 

We next evaluate the condition that <3?2 must be square integrable. Near x = we 
can use the relation 2-Pi(«> b, c, 0) = 1 to obtain the expansion 

$ 2 = (a 2^'+! + x i+i ) (1 + 0(x)) , 

showing that the parameter a must vanish. We conclude that $2 must be of the form 

$2 = /? (i-y) i+ Hv) i+ * 2 F 1 (-X + j + 1, X+ j + 1,2 + j,y), (A.19) 

To avoid trivialities, we can assume that the remaining parameter (3 is non-zero 
(namely, otherwise $2 would vanish identically and, according to (|A. 14|) . <3?i would 
also be identically zero). We then obtain a non-square-integrable pole at x = 71 unless 
the corresponding hypergeometric function 2 -Pi in ()A.19P vanishes at y = 1. From the 
expansion [TJ eq. 15.3.10] one sees that this is the case only if 



ro + i-A) ro + i + A) 



Since the T-function has no zeros (see [U eq. 6.1.12]), this condition can only be 
satisfied if the denominator is singular. Using that the Gamma function T(x) is only 
singular at the negative integers x = 0, —1, —2, . . . (see [TJ after eq. 6.1.3]), we obtain 
the following necessary condition. 

Lemma A. 4. The radial equations (lA.llj) admits a square-integrable solution $ £ 
L 2 ((0, 7r)) 2 only if X is of the form 

X = ±(n + j + 1) with n G N . (A.20) 

The result of this lemma greatly simplifies our formulas, because for these admissible 
values of A, the hypergeometric function in (|A.17p can be expressed in terms of a 
Jacobi polynomial (see [U eq. 15.4.6]), 

Mx) = (3{-^) {—^) ^TWn P " {cosx) 



(3n\ 



2i+l (2+j) 



(1 - u)2 + t (1 + u)2+i P^' +1 ^(n) , (A.21) 



where (2 + j) n denotes the Pochhammer symbol, and in the last line we set u = cosx. 
Substituting this result into (|A. 14|) we can compute the derivative using [1] eq. 22.8.1]. 
Furthermore, we can use [H eqns. 22.7.17-22.7.19] to modify the prefactors (1 ± u) 



20 F. FINSTER AND M. REINTJES 

as well as the upper indices of the Jacobi polynomials. By suitably applying these 
relations, one obtains 

*l(x) = ± 9 , + w 9 n i r. (1 - u)^ (1 + «)4+f P^ +l \u) , (A.22) 

where the signs ± refer to the two choices in ()A.20p . The L 2 -norm of the func- 
tions (IA.21j) and (jA.22j> can be computed with the rule (see [U eq. 22.2.1]) 



/_ 



'l-.r(i + ,)'(^(,^- 2 ° +M r(n + « + i,r(n + 6 + i) 



_i y v ' V " K ') 2n + a + b+l n\(n + a + b)\ 

In particular, one sees that <I>i and $2 are both square integrable. Choosing the 
constant such that the L 2 -norm of $ (as given by (|A.13p ) equals one, we obtain 

$A ^n!(n + 2i + l)! . . 1 / cos( X /2) Pi ij+1) (cos X ) 

1 > - v sm^+2( x ) . (A.23) 



±J> 2?+* r(j + 2) \± sin( X /2) Fi i+lj) (cos x ) 

We finally summarize our results. 

Theorem A. 5. T/ie intrinsic Dirac operator on S 3 has the spectrum (I2.13p . Working 
in the representation (|A.1|) . an orthonormal eigenvector basis ip^jk w ^ ^ €= No, J €= 
No + I and A; G {— j, — j + 1, . . . , j'} is obtained by the ansatz (|2.15j) wit/i $ = ($1, $2) 
according to (|A.23|) . T/ie corresponding eigenvalues are given by (|2.14|) . 

Proof. Noting that V S 3 is an elliptic differential operator on the compact manifold S 13 , 
standard elliptic theory (see [11]) yields that V S 3 is essentially self-adjoint and has a 
purely discrete spectrum. Since in Lemma IA.3I we constructed a basis for the angular 
dependence, and in Lemma [A. 41 we determined all eigenvalues of the angular operator, 
the resulting eigenfunctions ipnjk are clearly a basis of L 2 (S' 3 ) 2 . Furthermore, according 
to Lemma IA.3I and the normalization after (|A.22j) , the basis vectors all have norm 
one. Thus it remains to show that the functions tp n j k are orthogonal. This follows 
immediately from the orthogonality of the angular functions x k < 1 ( see Lemma IA.3|) 

as well as the fact that eigenvectors corresponding to different eigenvalues of V s s are 
orthogonal. □ 

We finally derive a useful formula for the spectral projectors. For any eigenvalue A, 
we denote the projector on the corresponding eigenspace by E\. It can be expressed 
as an integral operator 

(E x ^)(x)= [ E x (x,y)ij( y )d f i y , (A.24) 
Js* 

where x,y £ S 3 and dfi is the volume form on S 3 . In terms of the orthonormal 

^njk 

S 3 



basis ipnik OI " Theorem IA.51 the kernels become 

W-l 

E ±w (x,y)= E < jk (^n jk (y) .... {A.25) 

n=0 k=-j 3 1 1 n 

where the bar denotes the adjoint spinor with respect to the scalar product defined 
pointwise on the spinors. This expression simplifies considerably if the second argu- 
ment is evaluated at the north pole, y = n. 
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Lemma A. 6. The integral kernels of the spectral projectors E\ in ()A.24|) satisfy for 
any n E No the relations 

£ ±(n+ 3)(z,n) = i n + 2) ' fooB(x/2) P n hl \x) Tic-*Mx/2) Pi U \ X )) , 
1 2 > 8tt2 r(n + |) V / 

where n denotes the north pole \ = 0, x G S 13 is parametrized by the coordinates (x, y); 
and cr x is defined in (|2.12p . 

Proof. In view of ()A.23j) and (|A.10p . V , njfc( n ) vanishes unless j = \. In this case, 
from (IA.23j) and [U eq. 22.2.1] we conclude that 



x^o smx V 7rn! x^o sinx 

Furthermore, in the case j = A, the 2-spinors (|A.7j) simplify to 

Xi i = ~H= n > Xi i 



2 2 
1 



47r V / ' a" a V 1 



2 1 2 



1 / cos i? \ -| 1 (e- % *wv& 

Xi , i - r. — I „w „;„ „q ' Xi , i 



and thus 



E xLlX*! i=^l. E Xi + 1 X^ 

^— 2 2 2~2 4-7T ^— ' 2 + 2 2 _ 3 47T 



(n + 2)! 



7T n! 



fe=±3 fc=± 2 

Using these formulas in (|A.10|) and (1A.25I) , we obtain 

and substituting (jA23j) gives the result. □ 

Acknowledgments: We would like to thank A. Grotz and the referees for their careful 
reading and helpful comments on the manuscript. 
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